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A POPULATION OF HOT, DUSTY ULTRA-LUMINOUS GALAXIES AT Z - 2 

S.C. Chapman, 1 Ian Smail, 2 A.W. Blain, 1 R.J. Ivison 3 ' 4 
Accepted for publication in the Astrophysical Journal 

ABSTRACT 

We report spectroscopic redshifts for 18 //Jy-radio galaxies at mean redshift of z = 2.2 that are faint 
at both submmillimeter (submm) and optical wavelengths. While the radio fluxes of these galaxies could 
indicate far-infrared (far-IR) luminosities comparable to high-redshift submillimeter-selected galaxies (^ 
10 12 Lq), none are detected in the submm. We propose that this new population of galaxies represents an 
extension of the high-redshift submm galaxy population, but with hotter characteristic dust temperatures 
that shift the peak of their far-IR emission to shorter wavelengths, reducing the submm flux below the 
sensitivity of current instruments. Therefore, surveys in the submm waveband may miss up to half of 
the most luminous, dusty galaxies at z ~ 2. Mid-infrared observations with Spitzer will be a powerful 
tool to test this hypothesis. 

Subject headings: cosmology: observations 
starburst 



galaxies: evolution — galaxies: formation — galaxies: 



1. INTRODUCTION 



The microJansky radio population (Si.4GHz > 30 /iJy) 
has been the key to pinpointing and studying the submm 
galaxies - SMGs (Ivison et al. 1998, 2002; Barger, Cowie 
& Richards 2000; Smail et al. 2000; Chapman et al. 2001, 
2002a, 2003a). Approximately 40% of /iJy radio sources 
with optical magnitudes R > 23.5 (the optically faint ra- 
dio galaxies: OFRGs) are detected at Ssso^m ^ 5 mJy with 
SCUBA, and conversely 65-70% of SMGs brighter than 
this flux limit have reliable radio identifications. A spec- 
troscopic survey of radio-identified SMGs has measured 
redshifts for 73 SMGs (Chapman et al. 2003b, 2004a - 
C04), allowing us to constrain their dust temperatures, 
luminosities, star-formation rates, evolution, clustering 
strength, and dynamical and gas masses (Chapman et al. 
2003b, 2004a; Blain et al. 2004a, 2004b; Neri et al. 2003; 
Greve et al. 2004; Smail et al. 2003, 2004; Swinbank et 
al. 2004). With these measurements for a representative 
sample of submm-bright galaxies we can now study the 
properties of obscured galaxies at z ~ 2 as easily as the 
typically less-obscured UV-selected population (Steidel et 
al. 2004). 

However, there remains the question of the nature of the 
^60% of the OFRGs which are not bright submm sources. 
Chapman et al. (2003a) attempted a comprehensive study 
of the whole /iJy radio population, but could not constrain 
the nature of those OFRGs without submm detections. 
This was mostly due to a lack of spectroscopic redshifts 
for these galaxies, which could thus either be interpreted 
as moderate-luminosity star forming galaxies at interme- 
diate redshifts, z ~ 0.5, or alternatively, bolometrically- 
luminous galaxies at similar redshifts to the SMG popu- 
lation, but which are not detectable in the submm wave- 
band. The latter possibility exists because submm flux is 
a relatively poor proxy for the bolometric luminosity of a 
galaxy, being strongly sensitive to the characteristic dust 



temperature (Blain 1999; Eales et al. 2000). The depen- 
dence for galaxies at z ~ 2 follows the approximate form 
S850pm — Tj ' Ltir - even a small increase in Td implies a 
large decrease in observed submm flux density (see Fig. 1). 
To investigate the nature of the submm-faint OFRGs, 
we have undertaken a spectroscopic survey of this popula- 
tion in parallel with the SMG redshift survey of C04. In §2 
we describe the sample and our observations. §3 presents 
our findings and §4 discusses these and gives our conclu- 
sions. All calculations assume a flat, ACDM cosmology 
with n A = 0.7 and H = 72kms~ 1 Mpc" 1 . 

2. SAMPLE AND OBSERVATIONS 

The sample studied here is defined by selecting all ra- 
dio sources in seven separate fields (CFRS-03, Lockman- 
Hole, HDF, SSA13, Westphal-14, Elais-N2, and SSA22) 
which have either been targeted with SCUBA (Holland 
et al. 1999) in photometry- mode or lie within existing 
SCUBA maps in these regions and which lack an op- 
tical counterpart brighter than R — 23.5. The latter 
criteria has the effect of eliminating the lowest redshift, 
low-luminosity sources as well as optically bright AGN at 
high redshifts (the R > 23.5 condition implies that L* 
quasars and optical-AGN are excluded from our sample 
for z ~ 2 - Boyle et al. 2000). For comparison to earlier 
work on OFRGs we note that the typical color of galax- 
ies at R > 23.5 is (R-I) ~ 0.5 (Smail et al. 1995). We 
note that a less stringent cut in optical faintness (R > 23) 
leads to a significantly larger number of broad-line QSOs 
and low-redshift star-forming systems in the sample. 

The radio data in CFRS-03, Lockman-Hole, Westphal- 
14, Elais-N2, and SSA22 were obtained and reduced as 
part of our own programs (C04) and reaches r.m.s. sen- 
sitivities of 4-10 /zJy. Details of the Lockman-Hole and 
Elais-N2 radio data and their reduction can be found in 
Ivison et al. (2002), and this description is applicable to the 
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other datasets. The HDF and SSA13 radio data were ob- 
tained from E.A. Richards (private communication), and 
reach r.m.s. sensitivities of 8 and 5/iJy respectively. The 
HDF data is described in Richards (2000), while the SSA13 
data is so far unpublished (a subsequent reduction of the 
SSA13 data is described in Fomalont ct al. 2004). Submm 
fluxes were measured for all optically faint radio sources 
in these regions from either targeted photometry-mode ob- 
servations with SCUBA in these fields by our group (the 
technique is described in Chapman et al. 2001), or from 
our reduced SCUBA maps of these fields taken and re- 
duced from the JCMT archive, using a weighted 3-beam 
extraction. In the latter case to avoid contamination from 
bright submm sources unrelated to the OFRGs, the maps 
were first cleaned of all other sources brighter than 3<7 be- 
fore the fluxes were measured. Original presentations of 
some of the submm maps in these fields are given in the 
following papers: CFRS-03 (Webb et al. 2003), Westphal- 
14 (Eales et al. 2000), Lockman-Hole and Elais-N2 (Scott 
et al. 2002), HDF, SSA13 and SSA22 (Barger et al. 1999, 
2000). 

We then refine our sample further by discarding any 
sources which are either nominally detected in the submm 
at > 2.5cr or have 2.5-ct limits on their submm fluxes which 
are consistent with them being brighter than 5 mJy at 
850 /an. This guarantees that the final sample are much 
fainter in the submm than the C04 sample, which has a 
mean submni flux of Ssso^m = 6.6 mJy, with two-thirds 
brighter than 5 mJy. As we show later, the typical submm 
flux for the OFRGs in our survey is only ~ 0.5 mJy and 
so most of these galaxies lie far below the confusion limit 
of current blank-field submm surveys. 

These selection criteria result in a parent catalog of 60 
submm- and optically-faint radio galaxies. We note that 
there is a comparably sized sample of OFRGs which are 
not formally detected in the submm, but whose submm 
limits are still consistent with detections, with an average 
Ssso/jm = 2.1 mJy, which are discussed in Chapman et al. 
(2004a). Note that of the 169 radio sources whose redshifts 
were used in the analysis of Chapman et al. (2003c), none 
satisfy the criteria for inclusion in this catalog, as they all 
have R < 23.5. 

A random subset of 36 of the 60 submm-faint OFRG 
were spectroscopically observed with LRIS (Oke et al. 
1995; Steidel et al. 2004) on Keck in several observ- 
ing runs throughout 2002, 2003, and 2004, under gener- 
ally good conditions (<~ 0.8-1.0" seeing). All the spec- 
tra cover the observed wavelength range from 0.3 /xm out 
to as much as 0.8 /an (depending on the slit position in 
the mask, and the grating used on the red arm of LRIS: 
4001/mm or 6001/mm). Exposure times were 1.5-4.5 hr, 
split into 30-min integrations. Data reduction followed 
standard techniques using custom iraf scripts. One- 
dimensional spectra were extracted and compared with 
template spectra and emission-line catalogues to iden- 
tify redshifts. All identifications are based on multi- 
ple features, most prominently the Lya line, along with 
weaker stellar/interstellar/AGN features and/or contin- 
uum breaks (Fig. 2). 

We obtain secure redshifts for 18 submm-faint OFRGs 
in our sample giving a spectroscopic completeness for this 
sample of only 50%. While relatively low, this complete- 
ness level is sufficient to elucidate some of the basic prop- 



erties of this population. We note that all the OFRGs 
in our sample from the HDF and SSA13 fields (including 
both the sources presented in Table 1, and the sources 
we exclude from the present sample because their submm 
limits are consistent with marginal detections) were in- 
cluded in the OFRG samples in Chapman et al. (2003a). 
However, spectroscopic identifications have not been pre- 
viously presented for any of these galaxies. The OFRGs 
in the remaining five fields are presented for the first time 
here. 

We stress that this secure spectroscopic sample is con- 
servatively restricted to those galaxies with the very best 
spectral identifications (with two or more reliable fea- 
tures). Of the remaining 18 nominally unidentified sources 
in our sample, many have either solitary bright emission 
lines, or several weak features which would put them into 
the same redshift range as those presented in this paper. 
The emission lines are likely to be either Lya, when no 
continuum was detected and the wavelength was <4500A, 
or [On]A3727 if continuum was detected on both sides of 
the line and the wavelength was ^4500A. We illustrate the 
spectra for the securely identified sample in Fig. 2 and list 
their observational properties in Table 1: their positions, 
redshifts, submm and radio fluxes and optical magnitudes, 
limits on their dust temperatures Td, total infrared lumi- 
nosities Ltir, and spectral class (divided into starburst 
(SB), AGN/SB or AGN). L T i R was calculated, based on 
the redshift and radio flux densities, K-correcting the syn- 
chrotron spectrum with an index a = 0.8 (Richards 2000), 
and assuming the local far-IR-radio correlation (Condon 
et al. 1991, Helou et al. 1985). The derived T d and L T ir 
are plotted in Fig. 3. 

Our spectral classes are derived from the UV spec- 
tral properties as follows: SB, no CivA1549 emission de- 
tected at > 3cr above the noise; AGN/SB, detectable 
CIVA1549 emission, but also robust interstellar absorption 
lines (most notably SIIA1264, OI/SIA1303, and CIIA1335) 
which would be heavily diluted/undetected if an AGN 
dominated the UV continuum emission; AGN, showing 
significant CivA1549 emission with no detectable interstel- 
lar absorption features. The Civ/Lya ratio has been used 
as an AGN diagnostic in previous studies of high-redshift 
galaxies. This ratio has a mean of 0.2 for the AGN sub- 
sample of z ~ 3 LBGs selected by Steidel et al. (2002), 
and a value of 0.12 for the composite radio galaxy spec- 
trum presented by McCarthy (1993). McCarthy (1993) 
and references therein argue that the ensemble properties 
of high-redshift radio galaxies (indications of a hard pho- 
toionizing spectrum, alignment of emission-line and radio 
major axes, and possible correlation between L ra( jio and 
LfOml) suggest they are predominantly AGN-powered. All 
of our OFRGs with AGN or AGN/SB classifications have 
Civ/Lya > 0.12. We also note that strong CivA1549 emis- 
sion is present in some classical luminous radio galaxies 
which also exhibiting interstellar absorption lines in the 
UV (e.g., 4C41.17 - Dey et al. 1997). Both of these facts 
mean that our classifications should be interpreted with 
caution. 

3. ANALYSIS AND RESULTS 

Our survey shows that the submm-faint OFRG are a 
high-redshift galaxy population, the mean redshift of our 
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sample is z — 2.2 with an interquartile range of ±0.3 and 
a full range spanning z = 0.9-3.4. At such high redshifts, 
these galaxies are bolometrically luminous, all but one 
with L TIR > 10 12 L Q , and 4 with L TIR > 10 13 L Q (Ta- 
ble 1, Fig. 3), assuming the local far-IR-radio correlation 
holds. These luminosities are comparable to those esti- 
mated for the submm-detected population uncovered by 
SCUBA (C03), as shown in Fig. 3. Moreover, the red- 
shift range populated by the submm-faint OFRG is the 
same as that inhabited by SMGs (C03; C04), which have 
a median of z — 2.3 ±0.4. The two populations also 
have very similar optical and radio characteristics: the 
18 spectroscopically-identified, submm-faint OFRGs have 
a median i?-band magnitude of R = 24.6 ± 0.2 and a me- 
dian radio flux of Si.4GHz = 79 ± 26 ^Jy, comparable to 
the values seen for the submm-detected radio sources in 
the C04 sample: R = 24.1 ± 0.2 and Si. 4G Hz = 74 ± 6/xJy 
(although with the caveat that SMGs are not restricted to 
i?>23.5). 

The similarity of the OFRG and SMG populations ex- 
tends to their UV spectral properties, where a compari- 
son of Fig. 2 and the spectra presented in C03/C04 shows 
that the spectral characteristics of both submm-bright and 
submm-faint radio galaxies are very similar. The spectra 
of the submm-faint OFRG span a range from pure star- 
burst (similar to the z ~ 2-3 population selected in the 
rest-frame UV - Shapley et al. 2003) to low luminosity nar- 
row line or Typc-II AGN with enhanced NvA1240 and/or 
CivA1549 emission. Although we see no broad-line Type-I 
AGN, this is not particularly surprising given our opti- 
cal limit was chosen to exclude them. The galaxy spectra 
in which AGN signatures dominate comprise about 20% 
of the submm-faint OFRG sample, with the SMG sam- 
ple of C04 exhibiting a similar fraction. These AGN are 
spectrally similar to the rcstframe UV-sclcctcd AGN at 
z ~ 3 described in Steidel et al. (2002), although those 
galaxies are typically undetected in the radio. A further 
quarter of our OFRGs exhibit interstellar absorption fea- 
tures in the UV-continuum, but also show signficant AGN 
emission lines (e.g. CivA1549, HciiA1640), which we clas- 
sify as hybrid AGN/SB. The UV spectra of the remain- 
der of our spectroscopically-identified OFRG, consisting of 
half of the sample, show no detectable signs of AGN. An 
Ha survey of OFRGs and SMGs (Swinbank et al. 2004) 
broadly supports the pure-starburst classification from the 
UV; galaxies identified as SB from the UV typically show 
narrow Ha linewidths (<500km/s) and low [Nil] /Ha ra- 
tios typically of star-forming galaxies. Moreover, we note 
that those galaxies for which we failed to identify a red- 
shift were either undetected on the slit, showed no obvi- 
ous features in the weakly detected continuum, or showed 
only a single spectral feature as described above. Some 
of these galaxies are likely to be high redshift starbursts, 
since AGN lines would likely have been identified even with 
the faint continuum magnitudes. However, a low redshift 
starburst interpretation would still be consistent with sev- 
eral of these galaxies. We can conclude that less than half 
of the OFRG sample are likely to show any spectral signa- 
tures of an AGN, similar to the spectroscopically-classificd 
SMG population. 

In fact, there is only one characteristic which differs 
substantially between our submm-faint OFRG sample and 
the SMGs in C03/C04 - the average 850-/xm flux of the 



galaxies. The C04 sample of SMGs has an average flux of 
6.6 mJy, whereas the submm-faint OFRGs have a variance- 
weighted average of only Ssso^m = 0.5 ± 0.3 mJy, an order 
of magnitude fainter. Why are these apparently ultra- 
luminous galaxies, which share many features with the 
submm-bright, high-redshift SMG population so faint in 
the submm waveband? As illustrated in Fig. 1 this could 
cither reflect: enhanced radio emission, relative to the far- 
infrared, compared to that normally seen in star-forming 
galaxies; or a hotter characteristic dust temperature than 
that of the submm-detected population at these redshifts. 
We discuss these alternatives in turn. 

Radio emission from an AGN is one possible route to in- 
crease the radio fluxes of these galaxies, but not increase 
their far-infrared luminosities - and thus leave them unde- 
tectable in the submm waveband. We see some signatures 
of an AGN in half of our sample with robust redshifts, 
however to significantly perturb the radio fluxes of these 
systems we would require a substantial contribution from 
the AGN and hence a luminous, central AGN. The spec- 
tral features we see can arise at relatively large radius (for 
instance in the outer parts of the accretion disk around the 
AGN and beyond - e.g. Hutchings et al. 1998) and so it 
is possible that dust in the very central regions is obscur- 
ing the true luminosity of the AGN in the restframe UV. 
Optically-thin radio emission could still escape from these 
regions and hence high spatial resolution radio observa- 
tions would identify a strong central point source. More- 
over, X-ray emission from the AGN can also escape, and 
hence sensitive X-ray observations can be used to search 
for luminous, but highly-obscured AGN in these galaxies. 
Both of these observational tests can be applied to the 
OFRGs lying in the well-studied HDF region. 

Of the four OFRG lying in the HDF region, two have 
restframe UV spectra which are classed as starbursts, one 
is an AGN/SB and one is classed as an AGN. None of 
the galaxies spectrally-classified as starbursts or AGN/SB 
are detected in the 2-Ms Chandra HDF observation, while 
the AGN is (Alexander et al. 2003). To provide a more 
general comparison we take advantage of the spectral sim- 
ilarity of SMG and submm-faint OFRG to combine our 
sample with that of C04 to provide 12 luminous, dusty 
high-redshift galaxies with AGN signatures and 20 with 
starburst spectra within the Chandra HDF. In the 2-Ms 
Chandra image, the 12 sources showing AGN signatures 
in their UV spectra are all detected in the X-ray data, 
whereas only 9 of the sources showing starburst spectra 
are detected. Moreover, the average 2-8 keV X-ray flux 
of those SMGs/OFRG with AGN-like spectra are roughly 
an order of magnitude greater than the galaxies with star- 
burst spectra, whose X-ray fluxes are more consistent with 
that expected from the X-ray binary star emission from a 
~ 10 3 Moyr^ 1 starburst (see Alexander et al. 2004). 

The high-resolution radio observations of the sources in 
the HDF come from the MERLIN/VLA map of this re- 
gion (Muxlow et al. 2004; see also Chapman et al. 2004b) 
and have 0.3" resolution - sufficient to identify nuclear ra- 
dio sources on sub-kpc scales. These data show that the 
radio morphologies for the starburst-classified OFRG are 
often extended and clumpy on ~ 1" (~ 8kpc) scales, while 
some galaxies with signs of AGN in their UV spectra also 
show similar radio morphologies (Chapman et al. 2004b). 
This indicates that the majority of the radio emission from 



these galaxies is not coming from a central point-source or 
from low-luminosity analogs of classical radio jets. Hence, 
we conclude that the majority of these galaxies do not 
have AGN-enhanced radio emission and thus they must 
be bolomctrically luminous systems. 

Why then arc these galaxies undetected in the submm, 
where we detect comparably luminous dusty galaxies at 
these redshifts (C04)? The answer must lie in the charac- 
teristic temperature of the dust emission in these galax- 
ies. This is shown in Fig. 3, which depicts the distri- 
bution on the T(J-Ltir plane of the local IRAS galaxy 
distribution (Chapman et al. 2003c) and the locations of 
the radio-detected SMGs at z > 1 from C04. A repre- 
sentative flux limit for current submm surveys is shown 
by the shaded region. This demonstrates that hot, high- 
luminosity dusty galaxies could be at similar redshifts to 
the SMGs, z > 2, and detectable in the radio, but not 
at 850 ^m with SCUBA (see also Blain ct al. 2004b). Us- 
ing the average, error-weighted flux for the submm-faint 
OFRG sample we can estimate the characteristic dust tem- 
perature these galaxies must have: J£ 50 K, rather than the 
~ 36 K of the submm-dctectcd population. While the dust 
temperatures distributions of the SMG and OFRG popu- 
lations overlap, Fig. 3 shows that the samples are distinct 
in Td/(l+z). The total range in Td spanned by luminous 
galaxies at high redshift could be a factor of two larger 
than that implied by submm-selected galaxies alone (Blain 
et al. 2004a, C04). 

We conclude that the submm-faint OFRGs are most 
likely hotter relatives of the SMGs detected at high red- 
shift. These galaxies have far-infrared luminosities of 
>10 12 Lq, and appear to be predominantly powered by 
starbursts, as shown by the relative weakness of AGN sig- 
natures in the UV, X-ray or radio wavebands. While our 
Ltir estimates rely on the validity of the local far-IR-radio 
correlation (Condon et al. 1991, Helou et al. 1985), if the 
correlation is very different at z ~ 2, then it would affect 
the luminosity calculations for the SMGs and submm-faint 
OFRGs in a similar manner, and they would still remain 
comparably luminous populations. 

The volume density of submm-faint OFRGs at z ~ 2.2 
are comparable to that of the SMGs (C03, C04), with 
p = (6.2 ± 2.3) x 10~ 6 Mpc -3 for 1 < z < 3. This sug- 
gests that around half of the ultraluminous infrared galax- 
ies at z ~ 2.2 are missed by existing submm selection cri- 
terion, and that the total volume density of > 4 x 10 12 L Q 
galaxies at this epoch is 1.3 x 10~ 5 Mpc -3 . We illus- 
trate this by calculating the star formation rate densities 
(SFRD) for our sample of 18 OFRG, divided into three 
redshift bins (Fig. 4). We note that the far-IR emission 
from between 20-50% of our sample may include a con- 
tribution from an AGN, however, this down-ward revi- 
sion in the estimated SFRD would be compensated by 
including those (predominantly star- forming) OFRG's for 
which we have only obtained single-feature redshifts and 
which are therefore not included in this estimate. We 
translate our luminosity-density measurements into a star- 
formation density using the average standard calibration 
of 1.9xlO 9 L (M yr" 1 )- 1 (Kennicutt 1998) and add the 
results to the SFRDs of submm galaxies in Fig. 4. This 
reveals that the observed OFRG and submm galaxy sam- 
ples have a SFRD at z ~ 2 comparable to rcstframe-UV 
selected galaxies after correcting the latter for dust ex- 



tinction by a factor of five. The space density of these 
luminous galaxies has decreased a thousand-fold over the 
last llGyrs, to just <~ 10~ 8 Mpc~ 3 at the present-day 
where the IRAS survey selected all the most luminous, 
low-redshift, dusty galaxies regardless of dust tempera- 
ture. This evolution is far stronger than is seen in less 
luminous galaxies and underlines the importance of trac- 
ing obscured activity to understand the earliest phases of 
galaxy formation in the Universe. 

We can further demonstrate the decoupling of the un- 
obscured and obscured surveys of distant, star- forming 
galaxies by investigating the limited extent to which the 
huge infrared luminosities of the OFRGs are traced by 
their rest-frame UV emission. By design these galaxies 
are faint in the optical, and typically very faint in the 
UV. Using the deep UBR photometry of OFRGs lying in 
the HDF, Lockman, and SA22 fields, we determine that 
there is some overlap of the colors and UV luminosities of 
the OFRGs with populations of z ~ 2 galaxies selected in 
the UV (Steidel et al. 2004). Hence, some OFRGs have 
similar UV-predicted SFRs to the Steidel et al. galaxies 
(Reddy & Steidel 2004). However, these are well below 
the star formation rates implied by the radio emission 
(~ lO 3 M0yr _1 ). A similar mismatch arises locally when 
one considers the faint UV emission and blue colors of 
some dusty ULIRGs (e.g. Goldader et al. 2002). 

4. DISCUSSION AND CONCLUSIONS 

Looking at the similarity in the median redshifts of the 
SMG and submm-faint OFRG samples it is tempting to 
conclude that the evolution of these two populations must 
be similar, peaking strongly in SFRD at z ~ 2.3. How- 
ever, there are substantial differences in the K-corrections 
between submm- and radio-selected samples. While the 
SMG redshift survey (C03/C04) includes a radio-selection 
criteria, the total SMG population is constrained by both a 
knowledge of the underlying parent population (C04) and 
the benefits of a negative submm-wave K-correction that 
maps a flux limited survey into a luminosity-limited sur- 
vey. Unfortunately, at z £, 3 the radio fluxes of the submm- 
faint OFRG begin to fall below the current sensitivity lim- 
its of the VLA. Hence for the submm-faint OFRGs we are 
unable to state whether the population seen at z ~ 2-3 
can extend to much higher redshifts with a comparable 
luminosity density (analogous to the apparently constant 
comoving luminosity density in UV-selected populations 
at z > 2-4, Steidel et al. 1999; Giavalisco et al. 2004). 

Understanding the relationship between the submm- 
faint OFRGs, the submm-brighter galaxies (C04), and 
other high- redshift populations (Franx et al. 2003; Daddi 
et al. 2003; Steidel et al. 2004) will be helped by comparing 
the strength of the clustering, both within and between the 
different classes. Unfortunately, this will require substan- 
tially larger samples than are currently available. How- 
ever, it is worth pointing out that the statistics of pairs 
of submm-faint OFRG within 1200 km s -1 of each other 
in these fields are comparable to those found for SMGs 
by Blain et al. (2004a), and similar calculations for the 
OFRG sample would yield an estimate of the correlation 
length ro ~ 8 h[" o Mpc, amongst the highest clustering 
amplitudes seen for any high-z population. In addition, 
several of these submm-faint OFRGs lie within the SMG 
associations described in Blain ct al. (2004a), suggesting 
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a close relationship between SMGs and OFRGs. These 
arguments provide circumstantial support for our claim 
that we have indeed found a new population of extremely 
luminous galaxies, comparable to the SMGs. 

Our discovery of this new, luminous galaxy population 
at z ~ 2 could have several implications. The most impor- 
tant is that the census of the most luminous star formation 
at the peak epoch in quasar activity (z ~ 2.2) is seriously 
incomplete when only the submm-bright fraction is consid- 
ered. We suggest therefore that the z ~ 2 star formation 
rate density has been underestimated in studies to date 
(Steidel et al. 1999; Smail et al. 2002; C03). At least out to 
z ~ 2, the radio selection provides us with a less biased as- 
sessment of the total energetic budget. This conclusion is 
not changed qualitatively if we remove the half of the new 
sample that exhibit signs of an AGN in their UV spectra 
(especially if we consider that the half of the parent sample 
which have only single- feature redshifts, and hence are not 
included here, are likely to be starburst systems at similar 
redshifts). However, the current sensitivities achievable in 
the radio cannot detect such hot, luminous galaxies be- 
yond z ~3, and we are therefore not in a position to assess 
the incompleteness of current estimates of the obscured 
star formation rate density at higher redshifts. Observa- 



tions from Spitzer and the extended- VLA may complete 
the picture initiated in this study. 
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Fig. 1. — The far-IR/radio spectral energy distributions (SEDs) for galaxy templates with identical far-IR luminosities but characteristic 
dust temperatures corresponding to 30, 40 & 50 K greybodies with dust emissivities (3 = 1.60, 1.55 and 1.50 respectively (Dale et al. 2001). 
The points represent a typical 6 mjy SCUBA galaxy at 850 and 450^tm (dots) from the SMG sample of C04 with a dust temperature of 40 K. 
As can be seen, galaxies with similar luminosities and similar radio fluxes, but with SEDs which are characterised by a dust temperature 
hotter than ~ 50 K will have 850-/im fluxes less than ~ 1 mjy at z = 2 and thus will be undetectable in current submm surveys. Note that 
24/im emission is definitely detectable by Spitzer at a level of 0.1 mjy. 
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Fig. 2. — Spectra of all sources from our sample showing the typical spectral characteristics of our observations and this class of galaxies. 
The spectral features used to measure the redshifts and classify the spectra, either starburst (SB) or AGN, are marked. The identified lines 
are: OviA1038, Lya A1215, NvA1240, SioA1264, Cm]A1296, CnA1335, SiivA1396, CivA1549, HeiiA1640, and (for RGJ105204.22+572745.7) 
HcllA3203, [On]A3727, and Ca H&K. We use dashed lines for primary line IDs and dotted ones for supporting IDs. The reader should note 
the close similarity between these spectra and those measured for the SMG population by Chapman et al. (2003b, 2004a). 
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Fig. 2. — (cont) Spectra of all OFRG from our sample. 
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Table 1 
Properties of submm-quiet OFRGs 



• Td limits derived using 2<r limits from the submm. 

• Ltir was calculated assuming the local far-IR/radio correlation (Helou et al. 1985), with a TIR color correction term calculated 
at the Td limit. 
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Fig. 3. — The distribution of dust temperature, T^, versus total infrared luminosity (Ltir, 8/im-1100/^m) for submm-faint radio sources 
with spectroscopic redshifts from our sample. A lower limit on Tj has been calculated assuming a 2-c upper limit to the submm flux of each 
source. For comparison, we show the radio-submm galaxies from Chapman et al. (2003a, 2004a), with the average error bar shown to the 
right. The typical sensitivity limit of surveys precludes detection in the submm waveband of sources in the unshaded region. The Chapman 
et al. (2003b) derivation of the range of local IRAS galaxies from the 1.2-Jy 60/im catalog are shown as a ±2<r envelope. The average T<j 
for optically-bright radio sources at z = 0.3-1 from Chapman et al. (2003c) are shown as open squares, lying well within the local ±2cr 
distribution. 
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Fig. 4. — The evolution of the energy density (parametrized by SFRD) in the Universe with epoch. Submm measurements (Ssso > 5mjy) 
from Chapman et al. (2004a) are shown at the median value for each redshift bin. Radio-identified SMGs with spectroscopic redshifts are 
shown with solid squares, while an estimate for the 35% of SMGs undetected in the radio is shown with an open square. The new OFRG 
measurements are added to the submm-infercd fitted values to demonstrate the current lower limits on the extent of obscured star formation 
activity at high redshifts from populations that are fully quantified and detected. This observed radio+submm total SFRD is compared to a 
fit to the published estimates from optical/UV surveys corrected for dust extinction of the star formation density (short dashed lines show 
the drier envelope - the fit is derived by C04 from a compilation from Blain et al. 2002, and includes data from Giavalisco et al. 2004, Stcidel 
et al. 1999, Connolly et al. 1997, Yan et al. 1999, Flores et al. 1999, Yun, Reddy & Condon 2001). 



